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We investigate spatial and temporal cross-correlations between streamwise and normal 
velocity components in three shear flows: a low-dimensional model for vortex-streak inter- 
actions, direct numerical simulations for a nearly homogeneous shear flow and experimen- 
tal data for a turbulent boundary layer. A driving of streamwise streaks by streamwise 
vortices gives rise to a temporal asymmetry in the short time correlation. Close to the 
wall or the bounding surface in the free-slip situations, this asymmetry is identified. Fur- 
ther away from the boundaries the asymmetry becomes weaker and changes character, 
indicating the prevalence of other processes. The systematic variation of the asymme- 
try measure may be used as a complementary indicator to separate different layers in 
turbulent shear flows. The location of the extrema at different streamwise displacements 
can be used to read off the mean advection speed; it differs from the mean streamwise 
velocity because of asymmetries in the normal extension of the structures. 



1. Introduction 

Coherent structures are very effective in transporting momentum across velocity gradi- 
ents and thus contribute significantly to frictional drag in turbulent flows. Depending on 
the type of flow and the position of the layer being studied, different kinds of structures 
can be identified (Robinson, 1991; Panton 2001). In wall bounded shear flows, Robinson 
(1991) describes a dominance of streamwise vortices close to the walls and horseshoe-like 
structures in the outer region. For the intermediate region one might imagine a grad- 
ual transition in relative weight from one to the other. In transitional internal flows at 
low to intermediate Reynolds numbers streamwise vortices and streaks are also present 
(Eggels et al. 1994; Hof et al. 2004; Grossmann 2000), and a complete self-regenerating 
cycle for the dynamics, in which vortices drive streaks which then generate vortices 
through a shear instability, has been proposed (Waleffe 1997). The relation between 
internal and wall-bounded flow situations has been established through simulations in 
laterally confined geometries which show a similar dynamical behaviour (Hamilton et 
al. 1995). The presence of vortices and streaks in flows with homogeneous shear (Kida 
& Tanaka 1994; Schumacher & Eckhardt 2001) and calculations within rapid distortion 
theory (Nazarenko et al 2000) highlight the significance of the background shear for their 
evolution and dynamics. 

Streamwise streaks result from the mixing of fluid across the shear gradient as induced, 
for instance, by streamwise vortices. This is a linear process that suggests a causal relation 
between their appearance: the vortices have to be there first, and can then be followed 
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by the streaks. With the pointwise measurements in boundary layers in mind, we take 
velocity components as indicators for the structures: the streamwise turbulent velocity 
component u for the streaks and the wall- normal or shear component v for the vortices. 
The temporal correlation can then be verified for linear models of the vortex-streak 
interaction by direct calculation (Eckhardt & Pandit 2003). In particular, the model 
shows that the temporal cross-correlation function C vu (At) = (u(t + At)v(t)) t will be 
asymmetric, and it will have its extremum at a finite time-delay At. The vortex as 
measured by v has a chance to influence the streak in u for At > 0, but not for At < 0. 
The question we address here is the extent to which this causal relation is reflected in an 
appropriate temporal correlation function in fully developed turbulent flows. 

The correlations we are interested in can be obtained from two-point data, from mea- 
surements displaced in space or time or both. Apparently, the first such data were ob- 
tained by Blackwelder & Kovasznay (1972) in a turbulent boundary layer. Their correla- 
tions for y/6 ps 0.2 I show a weak asymmetry of under At — > —At. Later, Blackwelder & 
Eckclmann (1978), studied cross-correlations between normal derivatives of the velocity 
components, and confirmed the asymmetry, as well as a shift in the maximum towards a 
positive time shift. Lagrangian studies of this cross-correlation along particle paths also 
show this asymmetry (Pope 2002). 

Fully turbulent dynamics differs from the linear model with stochastic forcing by the 
presence of nonlinearitics and a self-consistent generation of the turbulent fluctuations. In 
addition, the self-sustaining cycle proposed by Waleffe (1997) for the complete dynamics 
consists not only of the non-normal amplification but also of a linear instability that 
generates normal vortices which arc then tilted into the streamwise direction by the 
main flow. The asymmetry is a property of the first part of the cycle, but not of the 
second. The correlation function of the full cycle will thus be a superposition of the 
contributions from the both parts, and the asymmetry will show up if the statistical 
weight of the second part is smaller than that of the first part. 

With this in mind, we want to turn to the analysis of the cross-correlation functions 
of the turbulent velocity components for evidence of the dynamical processes underlying 
the self-sustaining mechanism for the formation of coherent structures. We will study 
in detail data from a low-dimensional model, from direct numerical simulations (DNS) 
of shear flows and from measurements with hot-wire probes in a boundary layer. The 
model allows for a detailed tracking of the dynamics of the various contributions to the 
spatial and temporal correlations. The DNS allows for an extension to two-point cross- 
correlations in space and time since they do not have to rely on Taylor's frozen flow 
hypothesis. Finally, the experimental data, although restricted to time asymmetry, allow 
for much higher Reynolds numbers and for a systematic study of the dependence of the 
correlation functions on the distance from the wall. 

We take coordinates with x in streamwise, y in wall-normal and z in spanwisc direc- 
tions. The quantity we focus on is the correlation function between the normal velocity 
component v and the streamwise component u, displaced in the streamwise direction by 
Ace and in time by At, 



The averages are over time and also over all points in an x-z-plane at fixed height y (in 



C vu (At,Ax;y) 



(v(x, y, z, t)u{x + Ax, y,z,t + At)) x . Z:t 
{v(x,y,z,t)u{x,y,z,t)) x ^ tt 



(1.1) 



f Throughout this work, the boundary layer thickness S is defined as the distance between 
the wall and the height y where U = 0.99C/oo with the free-stream velocity [/oo. 
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Figure 1. The time correlation function for the four modal pairs that contribute to the 
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the model and the DNS). Time correlations at one point are given by 

(v(x, y, z, t)u{x, y,z,t + At)) x . z . t 



C vu (At;y) = C vu (At,0;y) 



(v(x,y,z,t)u{x,y,z,t))v 



(1.2) 



In order to quantify the asymmetry effects we introduce the following measure for the 
temporal cross-correlations 



Qvu(At) 



C vu (At) - C vu (-At) 



(1.3) 



C vu (At) + C vu {-At) 

(the dependence on height is suppressed in these expressions). For the extended correla- 
tions due to the non-normal amplification we expect ICV^Ai)! > \C VU {— At)\ such that 
Q vu > for these cases. 

The outline is as follows. In the next section, the low-dimensional model of a turbulent 
shear flow by Moehlis et al. (2004, 2005) is discussed. Section 3 describes the analysis of 
the DNS in a nearly homogeneous shear flow followed in section 4 by the discussion of 
boundary layer experiments of Knobloch & Fernholz (2004). Finally, a summary and an 
outlook are given. 



2. Low-dimensional model of turbulent shear flow 

The linear model of Eckhardt & Pandit (2003) can be extended to a nine-dimensional 
representation of shear flows, as discussed in more detail in Moehlis et al. (2004). The 
system is confined between two free-slip planes and driven by a volume force that sustains 
a laminar sinusoidal flow profile. The model captures the non-normal amplification pro- 
cess and completes it with modes for transversal shear and instabilities of the streamwise 
streaks. With L y = d/2, the aspect ratio is L x : L y : L z = 2ir : 1 : n, and we simulate 
the flow at a Reynolds number Re = U^dji^lv) = 180, where the reference value for the 
velocity Uq is determined from the sustained laminar velocity profile at y = d/4 with d 
being the distance between the free-slip planes. 

With the Galerkin modes u.;(x) and the amplitudes ai{t) we can write the turbulent 
velocity field as 

N 

u(x,t)=5^0i(t)ui(x). (2.1) 

i=i 

The spatial part of the cross-correlations can be calculated analytically from the pre- 
scribed modes u^(x). The temporal part follows from the numerical solution of a system 
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Figure 2. (a) Space-time cross-correlation C vu (At, Aa;; y) in the low-dimensional model taken 
at a height y — y/(d/2) = 0.52. The contour levels increase in steps of 0.1. The solid line 
represents the dimensionless mean velocity U(y = 0.52) = 0.74. (b) A cut through the contour 
plot in (a) at Aa; = gives the temporal cross-correlation C vu (At;y) . (c) Asymmetry measure 
Q vu (At) corresponding to diagram (b). This measure is independent of y due to the spatial 
mode dependencies. 

of ordinary differential equations for the dj(i) that results with (|2.1|) from the Navier- 
Stokes equations. Equation ll.ljl then becomes 

N 

C vu (At,Ax;y) = ^ ( ai (t) aj (t + At)) t ( Vi (x, y, z) Uj (x + Ax, y, z)) XjZ . (2.2) 

Of the 45 possible pairs of modes that could contribute to the correlation function only 
four terms, involving five modes, do. The modes that contribute arc: 112, a strcamwise 
streak with no streamwise variations; U3, a strcamwise vortex; U6 and 117, which describe 
two different wall-normal vortices; and Us, a mode that depends on all three coordinates. 
Thus C vu = C%u + C™ 6 + + where the superscripts indicate the particular mode 
couplings. Of these correlators, C^ 2 is exactly the one that probes the relation between 
streamwise vortices and streaks and therefore should have a significant variation with 
respect to At. This is indeed the Fig. shows. 

When all contributions are collected, the space-time contours for C vu (At, Ax; y) at 
y/(d/2) = 0.52 are obtained (Fig. . The time correlations at a fixed position correspond 
to a cut along Ax = (cf. Fig. lb) and the space correlations for fixed time from a cut at 
At = 0. The time correlation function is negative, as is to be expected for a flow where 
the streamwise velocity increases in the positive y-direction. It is asymmetric and the 
asymmetry gives a positive Q in the center of the layer. Because of the small number 
of modes in the system a more detailed dynamical system study is possible. An analysis 
of the periodic orbits in Moehlis et al. (2005) shows that some of them clearly follow 
the vortex-streak-instability sequence, but several do not. Nevertheless, the correlation 
functions in Fig. [2] show that the temporal asymmetry expected from the non-normal 
amplification process persists even after taking time averages. 

Besides this similarity, there are differences to the linear model. The contour levels 
(Fig. la) show a global minimum that is shifted from the origin toward (At, Ax) « 
(1,-1) for this particular height. The shift in time stems from the streak- vortex coupling 
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N x x N y x N z L x : L y : L z U S e u rms S* R\ 



DNS-1 128 x 65 x 128 2tt : 1 : 2tt 1 2 0.04 0.37 6.1 79 
DNS-2 256 x 129 x 256 2n : vr : 2tt 3/tt 6/tt 0.44 1.08 5.0 166 

Table 1. Parameters of the two DNS that were taken for the analysis. U is the mean 
streamwise velocity at the free-slip boundary. The mean energy dissipation rate is given by 
e = 15v{{du/dx) 2 ), as in the experimental determinations. S = dU/dy is the constant shear 
rate. The root mean square velocity is given by u rms = {u 2 } 1 ^ 2 , the dimensionless shear param- 
eter is S* = Su 2 ms /e and the Taylor microscale Reynolds number _Ra = \J 15/ (ve) u 2 ms . The 
spectral resolution criterion k max n > 1 is satisfied with k max = \> r 2N x /3. 



contributions C^ 2 and C^f, and is compatible with the stochastic model (Eckhardt & 
Pandit 2003). The one in position can be traced back to the coupling of mode 8 with 
itself, (vg(x, y, z)u%(x + Ax, y, z)) x ,z = 7r/(5 + n 2 ) sin(Ax/2) sin(iry) for — 1 < y < 1. The 
space-time contours are elongated along an axis whose slope has dimensions of velocity. 
This velocity is not the mean velocity at the height of the measurement. As we will show 
below, this is due to an asymmetry in the width of the structures in the normal direction. 

3. Nearly homogeneous shear flow 

The direct numerical simulations of a turbulent shear flow also refer to a flow bounded 
by two parallel free-slip plates, driven by a volume force that sustains a linear shear flow 
in the mean, U(y) = Sy, except for a small boundary layer near the plates. Details of 
the numerical simulations are given in Schumacher & Eckhardt (2000) and Schumacher 
(2004). Relevant parameters for the simulation are listed in table 1. 

Space-time contours for the cross-correlations in run DNS-1 are shown in Fig. 
Several features are similar to the ones in the low-dimensional model: it has the same 
asymmetry with respect to time and the iso-countours are oval and not aligned with 
the coordinate axis. Thus, even though more spatial degrees of freedom are present, 
the Reynolds number is higher, and the turbulence is fully developed, the non-normal 
amplification is reflected in the Eulerian cross-correlation function. 

However, there are noticable additional features. The asymmetry measure shows a 
pronounced height dependence, being strongest close to the bounding surfaces and getting 
weaker towards the center. In addition, it shows a time interval where its value is negative, 
Quv < (sec Fig. Et)- This interval is almost negligible close to the walls and becomes 
longer as the reference position moves towards the center. We see this phenomenon linked 
to the difference in the number of modes that can contribute, and hence to the possibility 
of additional dynamical processes. Assuming that the smallest scale is set by dissipation 
and does not vary much across the flow, the largest scale for the possible structures is 
set by the distance to the free-slip boundary. By this reasoning there arc fewer active 
modes close to the wall than in the center, thus limiting the nonlinear interactions and 
highlighting their correlations. 

A spatial plot of the streamwise turbulent velocity component u and the shear compo- 
nent v in the x-z plane for one point in time and at a fixed height y in the layer (see Fig.01 
reveals that the contributions to the cross-correlation function come from fragmented re- 
gions, of an extension compatible with the dimensions of coherent structures. Negative 
contours of v indicate streamwise vortices which generate the streamwise streaks (shown 
as gray-filled contours of u). The maxima of u and v contours are displaced slightly, in 
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Figure 3. Space-time cross-correlation C vu (At,Ax;y) of a nearly homogeneous shear flow. 
Data are taken from DNS-1 (see table 1) at a height y/L y — 0.11. The Reynolds number 
Re = ULy/v = 1800 where U is the mean turbulent velocity at the boundary, (a) Space-time 
plot of the cross-correlations. The contours increase in steps of 0.1 and the unit value is at the ori- 
gin. The solid line represents the dimensionless mean turbulent velocity U(y/L y = 0.11) = 0.77. 
(b) Temporal cross-correlation along Ax — for different heights: A : y/L y = 0.13, 
: y/Ly — 0.25, o : y/L y — 0.50. (c) Asymmetry coefficient for C vu (At;y) at the same 
heights as (b). 



accordance with the off-set in the maximum of the spatial cross-correlation in Figs. 2a 
and 3a. 

The off-set can be explained by the observation that a streamwise vortex pair centered 
at height y will be advected with the corresponding mean streamwise velocity at that 
height, U(y). The pair will mix slower moving fluid into a region that moves on average 
faster, and hence will temporarily reduce the local advection velocity U c to values below 
the mean velocity U(y). But since U c will advect the streamwise streak that is about to 
be lifted up, it will remain behind the vortex pair, resulting in the spatial shift of the 
most intense cross-correlation. 

The correlation functions shown in Fig. 2 and 3 and many others for different aspect ra- 
tios and Reynolds number show an inclination of the isocountours in the spatio-temporal 
correlations C vu (At, Ax; y). Since the two axes being compared have dimensions of time 
and length, the inclination has the dimension of a velocity: but as the comparison with 
the straight lines in both figures shows, the velocity with which these structures are ad- 
vected is systematically lower than the mean velocity at that height, U(y). We can trace 
this effect back to an asymmetry of the spatial autocorrelation function of the streamwise 
turbulent velocity in the normal direction, as measured by 

C uu (Ay; y ) = (u(x, y Q , z, t)u(x, y Q + Ay, z, t)) x ,z,t ■ (3.1) 

The left and middle diagram of Fig.[S]show that this function is asymmetric with respect 
to the wall-normal direction, obviously influenced by the presence of the free-slip walls 
at y = and y = L. If we estimate the mean advection speed of the coherent structures 
from an average of the streamwise speed over a domain determined by the full width at 
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Figure 4. Snapshot of the two turbulent velocity fields entering the cross-correlation function. 
A slice cut from a sample of DNS-2 (see table 1) is taken. The turbulent streamwise velocity u 
is indicated by shading for values of u between and 1.5, only. The contour lines show three 
isolevels of the wall-normal component v at values of —0.8, —0.5 and —0.2. The maxima are 
shifted relative to each other by a small streamwise distance that corresponds with the shift of 
the maximum of the space-time cross-correlation by Aa; as visible in Fig. 3a. 

half maximum of C uu (Ax, Ay; yo), wc find 



Here, £2 and £\ are the widths at half maximum of the asymmetric C uu (Ax = 0, Ay; yo) 
(see also the mid diagram of Fig. |3J) . The result of such an averaging procedure can be 
seen in the right diagram of Fig. |SJ The convection velocity as defined by l|3.2[l becomes 
smaller as the mean velocity and coincides with the inclination of the space-time contours 
of the velocity cross-correlations of Fig. 3. The coherent structures thus move with the 
streamwise speed as determined by an average over their size. 

4. Turbulent boundary layer 

The third class of flows for which we determine cross-correlation functions are high- 
Reynolds number boundary layer flows. The experiments were done at a wind tunnel 
of the Hermann-Fottingcr Institute (HFI) in Berlin and the German-Dutch Windtunnel 
(DNW). Triple hot-wire probes allowed measurements of all three velocity components. 
With sampling rates of 20 kHz at HFI and 125 kHz at DNW, data sets containing 
about a million data points at the HFI and 5 million at the DNW could be obtained. 
The boundary layer thickness (see the footnote on p. 3) was found to be i5 = 63mm 
at HFI and S = 240mm at DNW. Some parameters are summarized in table 2; further 
experimental details may be found in Knobloch & Fernholz (2004). 

Because of the measurements at a single location, only the short time behaviour of the 
cross-correlation functions can be determined. The results for different distances from 
the wall and different Reynolds numbers are shown in Fig. El Already for the data set 
closest to the wall there is a time interval with negative values in the asymmetry measure, 




(3.2) 
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Figure 5. Left: Contour plot of spatial autocorrelation function C uu {Ax, Ay; yo) as given by 
l|3.1|l in streamwise and wall-normal directions taken at j/o = 1/4 from DNS-1. The inclination of 
this structure is about 9 degrees. Mid: A cut through the contour at Ax = is shown. It is also 
indicated how the two widths at half of the maximum, l\ and £2, of the resulting autocorrelation 
function are defined. Right: Comparison of the mean (U) and the convection (U c ) velocities for 
data from DNS-1. Black solid line is for mean velocity as a function of the wall-normal distance 
y. The symbols stand for the corresponding convection velocity as given by 13.21 . 



and this interval increases as one moves further out. Comparison with the homogeneous 
shear flow DNS suggests that this point is already in the transition region away from the 
vortex dominated near wall layer. For the points furthest from the wall no reversal to 
positive values is detected. 

The data from HFI and DNW are collected at about the same relative positions when 
heights are measured in units of the boundary layer thickness, y/S. The asymmetries 
at these heights show remarkably similar behaviour, especially for the value y = y / 5 = 
0.11, which is present in both data sets, and for the two similar values y = y/S = 
0.31 (HFI) and 0.34 (DNW). In wall units, the y = 0.11 corresponds to a height of 
y + ^ 174 for the HFI data and y + ~ 3 950 for the DNW data In a greement with 
the findings of be Graaff fc Eaton 1 l|2000h and iDel Alamo et al. I <|2003l) for turbulent 
intensities, the cross-correlations collapse in external scaling, i.e. relative to boundary 
layer thickness and the external velocity. A possible explanation for this behaviour could 
be that the intermittent bursting activity in the boundary layer lifts fragments of the 
coherent structures higher into the intermediate layer, where their further breakup is 
determined by the boundary layer thickness S (see e.g. Blackwclder & Kovaszany 1972). 
Analysis of additional data shows that the asymmetry in the time correlations shows up 
for positions up to y = 0.05<5. 

Information about the instantaneous in-plane correlations similar to Fig. [SJi can be 
obtained from PIV measurements at fixed heights. The data in Fig. [7] show that the 
inclination in C uu (Ax, Ay; y ) is preserved and has about the same value. 



5. Summary 

A comparison between the three sets of data shows that the vortex-streak interaction is 
reflected most strongly in the cross-correlation function closest to the wall or to bounding 
surfaces. Further away the signal gets weaker, as expected by the change in structures 
(Robinson 1991). The similarity of the asymmetry measure further out suggests the 
prevalence of a similar dynamics. Interestingly, the asymmetry measures are similar when 
the height is measured in units of boundary layer thickness, rather than viscous length 
scales. It will be interesting to see how horsehoe vortices and their dynamics are reflected 
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Table 2. List of boundary layer measurements. We have picked three distances from the wall 
for every free-stream velocity Uaa. y + = yu T /v with u T = (r wa ii / ' pf ■ The other quantities are 
defined as in the caption of Table 1. 
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Figure 6. Asymmetry coefficient Q vu (At) for two sets of turbulent boundary layer data. Time 
is given in units of S/Uoo for both figures, (a): HFI measurement at Reg = 41600 for three 
different heights y = y/5 above the wall, k : y = 0.02, ■ : y = 0.11 and o : y = 0.31. (b): 
DNW measurement at Reynolds number Reg = 1237900. Here k : y = 0.02, ■ : y = 0.11 and 
o : y = 0.34 (see also table 2 for more details). 



in correlation functions, and whether they can explain the cross-correlation functions or 
whether other dynamical processes have to be identified. 
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